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^ . Abstract. The imposition of symmetries or special geometric properties on sub- 

C^ I manifolds is less restrictive than to impose them in the full space-time. Starting from 

V/^ ' this idea, in this paper we study irrotational dust cosmological models in which the 

Cn I geometry of the hypersurfaces generated by the fluid velocity is flat, which supposes 

^^ ' a relaxation of the restrictions imposed by the Cosmological Principle. The method 

^ . of study combines covariant and tetrad methods that exploits the geometrical and 

^^ I physical properties of these models. This procedure will allow us to determine all the 

■^y ' space-times within this class as well as to study their properties. Some important 

Qs , consequences and applications of this study are also discussed. 
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:• ■ 1. Introduction 

^ • 

rS I One of the most interesting issues in modern cosmology is the determination of the 

P^. ', geometry of our (local) universe, an important piece of information to understand 

cosmological phenomena and to interpret correctly observational data from them. This 
question has been addressed in the literature many times and from many different points 
of view. On most occasions the starting point was the Cosmological Principle, which 
assumes homogeneity and isotropy and leads to the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) space-times, which admit a six-parameter group of isometrics whose 
surfaces of transitivity are spacelike hypersurfaces of constant curvature (see, e.g., [|l|). 
Hence, the question reduces to know whether the universe corresponds to the case of zero 
(fiat models), or positive (closed models), or negative (open models) spatial curvature. 
However, we know that the Cosmological Principle is an idealization and that some 
important cosmological issues, as for instance the description of the formation of cosmic 
structures, require to consider the presence of inhomogeneities in the distribution of 
matter and therefore, they must be studied using inhomogeneous cosmological models 
(see Pi for a review). 

The main idea on which this paper is based is that we can relax the implications 
of the Cosmological Principle in the following sense: Instead of considering that the 
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restrictions coming from the homogeneity and isotropy assumptions are placed on 
the space-time manifold, we will consider such restrictions on certain submanifolds. 
More precisely, we will assume that the space-time is described by an irrotational dust 
solution of Einstein's equations, which seems to be a good approximation to describe 
the universe in the matter-dominated era, and then we will consider that the spacelike 
hypersurfaces orthogonal to the fluid velocity admit a six-parameter group of motions, 
i.e., a maximal group of isometries. The important point is that these symmetries are not 
necessarily symmetries of the full space-time, and therefore they are less restrictive. This 
kind of symmetries, those imposed on submanifold, have been usually called intrinsic 
symmetries P]. In cosmology, intrinsic symmetries have been studied by several authors. 
In particular, Collins and Szafron Q] studied them systematically in inhomogeneous 
cosmological models. Another interesting approach is the general programme put 
forward by Stephani and Wolf P, who made a systematic study of vacuum and perfect- 
fluid space-times admitting fiat hypersurfaces with special properties of their extrinsic 
curvature. 

As we have said before, in this paper we will study Irrotational Dust Models (IDMs 
hereafter) in which the hypersurfaces generated by the fluid velocity admit a maximal 
group of isometries, which means they have constant curvature. For simplicity, we will 
only consider here the case in which this group is such that the hypersurfaces are flat 
(vanishing curvature). The interest of this study is emphasized by recent observations P 
that point to the fact that the spatial geometry of the universe is close to be flat. On 
the other hand, this kind of models were already considered by Collins and Szafron [^ 
(Paper III). However, in their study they neither find exact solutions nor finish the study 
of the propagation of the constraints. In the present work, we have used a different point 
of view and a different procedure to study these models. Briefly, the plan of this paper 
is the following: In section ^, we will study these models from a completely covariant 
point of view by using the fluid covariant approach. In section |, we use the results 
of the covariant study to choose a preferred orthonormal basis, which will allow us to 
go further into this study and to find all the possible IDMs with flat spatial geometry. 
In section ^, we discuss the results and their main consequences, as well as some other 
related questions. 

2. Covariant approach 

Irrotational dust solutions of Einstein's equations are models suitable to describe the 
formation of large-scale cosmic structures as well as the dynamics of the Universe in 
the matter-dominated era. The energy-momentum distribution is described by the fluid 
velocity u and the energy density associated with it, p, the energy- momentum tensor 
being then 

Tab = pUaUb, U"-Ua = -1, p>0. (l) 

The fluid velocity is irrotational {u\^adhUc] = 0), i.e., it generates orthogonal spacelike 
hypersurfaces, and taking into account the conservation equations for the energy- 
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momentum tensor (|1|), we deduce that it is also geodesic {u^VhU"" = 0). Hence, there is 
locally a function r(x") which is at the same time the proper time of the matter and 
the label of the hypersurfaces orthogonal to the fluid velocity, which will be denoted by 
S(r). In terms of r we can write u as follows 

u = d-r , u = —dr . (2) 

And then, using adapted coordinates {r, y'^}, where y'^ are comoving geodesic normal 
coordinates {u'^day^{x^) = 0), the line element can be written in the form 

ds^ = -dr^ + h^,{r, y^)dy^dy'' . (3) 

In this paper we will study IDMs in which the hypersurfaces S(r) are flat. They are 
covariantly characterized by the vanishing of the Riemann tensor of these 3-dimensional 
manifolds 

'R\cd = 0. (4) 

This means that for each r we can find comoving coordinates y'^ such that the induced 
metric on E(r) becomes the fiat metric, i.e., /i/xi/|e{t) = S^lu■ 

Taking into account the properties of these models, the fluid covariant approach 
to relativistic cosmology is an adequate tool to study them. The variables used 
in this formalism are associated with the fluid velocity u, and in our case they are: 
The orthogonal projector to u, hab = gab + UaUb, which at the same time is the 
first fundamental form of the hypersurfaces S(r); the non-zero kinematical quantities, 
the expansion O = VaW" and the shear tensor aab = V(a'Ufc)||, which determine the 
second fundamental form of the hypersurfaces S(t), Kab = ^Qhab + o-ab] the energy 
density, p; and the gravito-electric and gravito-magnetic tensors, Eab = CacbdU^W^ and 
Hab = *Cacbdu''u'^ {E{ab) = Eab, H^ab) = Hab) respectively, where Cabcd denotes the Weyl 
tensor and * the dual operation. 

The equations governing these quantities can be obtained by sphtting, with respect 
to the fluid velocity, the Ricci identities applied to u, (VcV^ — VdVc)u"' = R'^bcdU^, 
and the contracted second Bianchi identities (equivalent in four dimensions to the 
non-contracted ones), V^Cabdc = "^[aTbic + \gc[a^b]T'^d (see, e.g., [0]). In order 
to write explicitly these equations we also need to split the space-time covariant 
derivative, Va, into a time derivative along the fluid worldlines, A°'"'b:- = u'^^cA"'"b--, 
and a spatial covariant derivative tangent to the hypersurfaces S(r), DcA"-"b-- = 
h°'e ■ ■ ■ h^ bhc'^V dA^" f-1 where A°'"'b... is any arbitrary tensor field. Moreover, it is very 
useful to introduce the spatial divergence and curl of an arbitrary 2-index tensor|| Aab 
(see [§ for details on the notation) 

div {A)a = D^Aab , curl Aab = ScdiaD'^Ab)'^ , 

X The angular brackets on indices denote the spatially projected, symmetric and trace-free part: 

A{ab) = [^(a'^^6) " 3'*'^ hab]Acd- 

§ These definitions are analogous to those for vector fields: div(y4) = DaA°- and curM^ = eabcD^A"^ . 
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where Eabc = Vabcd.u'^ {^"'''^^def = '^^■h^"'dh^eh^^ f) is the volume 3-form of the hypersurfaces 
S(r), and rjabcd the space-time volume 4- form. Finally, for two arbitrary spatial 
symmetric tensors, Aab and Bab, we define the commutator as 

[A, B]ab = 2A[a'^Bf,]c , [A, B]a = ^eabc[A, B] '^ = SabcA dB"^ . 

Using all these definitions, the resulting equations can be divided into two 
differentiated groups of equations. The first one provides the evolution of our variables 
along the fluid worldlines, and usually they are called the evolution equations. In the 
case of IDMs they look as follows (see, e.g., [§) 

hab = 0, (5) 

e = -i02-a"V„,-ip, (6) 

C^afe = — 3 0O"a6 " O'cia'^b)'' " -E'afe , (7) 

p = ~ep, (8) 

Eab = -QEab + ^(yc{aEbf - \pOab + CUrl Hab , (9) 

Hab = -QHab + 3(r,(^aHbf - curl Eab . (10) 

The other group of equations are relations between spatial derivatives of the variables, 
which are usually called constraint equations. In the case of IDMs they are (see, e.g., 

i) 

C\ = dwia)a-lDae = 0, (11) 

C'^ab = curl aab - Hab = 0, (12) 

C\ = diY{E)a-[a,H]a-lDaP = 0, (13) 

C\^dw{H)a + [a,E]a = 0, (14) 

C^ = p + |a"V,b - |e2 - i3i? = , (15) 

C ab = Eab — Sab + Crc<aCb>'^ — gQ^-afe = , (16) 

where ^R and ^Sab are the trace and trace-free parts of the Ricci tensor, ^Rab, of the 
hypersurfaces S(r). Here, it is worth to note that the last two constraints (|T5|,|l6|) come 
from the Gauss equations 

'R\,a = h\R-jg,hf,h\h\ - K\Km + K'^dKbc , (17) 

which relate the Riemann tensor of the hypersurfaces S(r), ^R"'bcd, to that of the 
space-time. Moreover, we have taken into account that for 3-dimensional Riemannian 
manifolds the Riemann and Ricci tensors contain equivalent information, and that the 
relationship between them is 

R°'bcd = — 2/l"be[c Rdf — 2^^bcd R {habcd = hacHd — hadhbc) ■ 
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Then, the consequences of (|) and (|T5|,|T6|) are 



P 
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1 ^ab-. 
2^ <^ab , 



's. 



ab 



=^ Eab = ^QcTab — Cy{a''(^b)c 



(18) 
(19) 



We can check that using these expressions for p and Eab the constraint (|1^) is 
automatically satisfied. Furthermore, from equation (|T^) it follows that the gravito- 
electric and shear tensors commute 



[<y,E], 



aacE' 



Ear-a" 



0. 



(20) 



And from the constraint ([1^) , we deduce that this is equivalent to the vanishing of the 
spatial divergence of the gravito- magnetic field (see, e.g., H, P, p!0(| ) 

div(i7)„ = 0. (21) 

m 



IDMs with a divergence-free gravito-magnetic field have been studied in 



In these 



works the consistency of the new constraints (pO| , ^lD under evolution was studied. The 
conclusion was that they lead to the following new constraint 

[a, curli/] = , (22) 

and that evolving repeatedly this constraint we would get an indefinite chain of new 
constraints. However, in the case of IDMs with fiat spatial curvature we have two 
additional constraints, namely ( p^p^Q]) . Using the evolution equations (|5|-p!0D it can be 
shown that the evolution of the first one, i.e. (|18[), does not lead to any new constraint, 
or in other words, it is consistent with evolution. On the other hand, it can be shown 
that the evolution of the second one, i.e. (^), leads to the following new constraint 

cnrlHab = . (23) 

In order to obtain this result we have used the following property of the shear tensor 



(7 (Tc<a<^b>d 



l{<y'"^crcd)<yab , 



which, in fact, is satisfied by any symmetric and trace- free tensor in 3 dimensions. As is 
clear, the new constraint (|23|) implies the constraint (p2[) . Therefore, the only constraint 



that we need to analyze in our case is (23). Before to do that, it is worth to note that 
I imply 



conditions (^ 



E)aHbc = D^aHi 



bc> 



which, following [|TT|, means that the spatial covariant derivative of the gravito-magnetic 
field has only the completely symmetric and trace-free part, which was called the 
distortion of Hat- In |Tl|, it was argued that this part must be non-zero if the space-time 
contains gravitational waves. 

Following our study, we have investigated the time evolution of the constraint (|23|). 
To that end, it is useful to introduce the following definitions 



Hnhr = D^nHl 



<aJ^bc> 



(^abc = D^aO'bO , Ta = D (7ba 



iDaQ. 
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and the decomposition in irreducible parts, with respect to the 3-diniensional rotation 



group, of the covariant derivative of the shear tensor |Tl 



Then, the result from the evolution of the constraint (P5|) can be written as follows 



+'^£cd<aHb>'^^0- e + 5H^a^Hb>c + —£cd<aHb>'^T + ^ {cTabcT'^ + (TabDcT^ — ac<aD'^Tby 

-\a<a'DbyT, - fj<aTb>} = . (24) 

As we can see, this expression cannot be reduced to an identity = by using the 
previous constraints, therefore it is a new constraint. Moreover, it only contains spatial 
derivatives of the shear tensor since T^ corresponds to its divergence and Hab [through 
the constraint C^ (12)] to its curl. The next step in the covariant analysis would be to 



study the subsequent evolution of this constraint. Taking into account that the evolution 
of the shear (^ can be written simply as [using (0)] 

C^ab = -Q(^ab , (25) 

and the form of the commutator between time and spatial derivatives 

[• , Da]A'-,... = -(|e5/ + a/)D,A'-,... + Hj {e,,'A'-,... + ■■■ + e,\A'\... + ■■■}, 

where A^'" c-.- is an arbitrary tensor, we deduce that all the relations we would obtain 
from ( P^ ) would be, like the constraints ( [23|J2^ ), of the form F [6, Oab-, DcCTab, DcDd<7ab] = 
0. However, to find these relationships would suppose a very big amount of calculations. 
Instead of this, in the next section we will follow the study using a preferred basis. 

3. Study in the preferred basis 

The important point now is the fact that from equations (p!9[) and (p5[) we can deduce 
that there is an orthonormal basis adapted to the fluid velocity u, {u, Ba} (n ■ Bq, = 0, 
^a-e/j = 6af3), such that the shear and gravito-electric tensors diagonalize simultaneously 

(T„/3 = Eaf3 = for a^ P, 

and all the basis vectors are parallelly propagated along u 

ii'^ = u^Vbu" = , e«" = u^Vbea'' = . 

This result was already shown in |^ for IDMs with a divergence-free gravito-magnetic 
tensor, i.e., satisfying (plj ). 

As is clear, the shear tensor has only two independent components in such a 
basis, which can be described by the following two quantities: o"+ = — 3o"ii/2 and 
o-_ = v^(cr22 — 0-33)72. From (|6ip^ and (pSf), the evolution equations for the expansion 
B and these two components of the shear can be written as follows 

= -1(92 + ^2+^2), 

a_^_ = -Qa+ , (26) 
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a_ = -0CT_ , (27) 

where from now on, Q = u{Q) = drQ for any quantity Q. The important point here is 
the fact that this system of evolution equations is cfosed, contrary to what happens in the 
general case, where p and Eat appear [see equations (QJ^)]. Using adapted coordinates 
{r, y"} [see (|])], we can integrate these equations and express the result in the form 

© = -, cr+ = ^— t, o-_ = ^— , V = -opir - ot){t - ot') , (28) 
V V V A 

where the subscript o on a scalar means that it does not depend on the proper time r. 

In these expressions, oCr+, oO"-, and qT are arbitrary, qt' is given by 

/ 4o(T 2_-'-/2, 2\ 

or = oT 7= — , oO- = -(oO-+ + oO-_j, 

V 3 oP 3 ^ 

and oP can be chosen arbitrarily since, as we can see from the equation (^), v is defined 

up to a proportionality factor depending only on ?/". Then, we can partially fix v by 

choosing oP to be a constant (we could also choose the value). Moreover, introducing 



(p8|) into (|T8D we get 

op 
P = —^ 

V 

so oP is directly related to the energy density. 

On the other hand, the evolution equations for the triad {£„} are 







>-i 


3VW 


^^+y>-i ; 








e2^ 


= -1(8 


+ 0-+ + V3cr_ 


-)e2^ 






es^ 


= 4(8 


+ cr+ — -\/3cr- 


-)e3^ 


which. 


using 


m, 


can be 


3olved to give 










= oe„'^(r-„r)-^"(r- 


-oT)-"", 



(29) 



where q^ = 2/3 — p^ and 

V3 oO" - 2 oCr+ o(T+ + V3(oCr + o(y^) oO'+ + V3(oO" - oO" 

Pi = ^^S^ ' P2 = —J= , P3 ~ 



3\/3oC 3v3oCr 3v 3 oCT 

We can check that Pa and go, satisfy the following relationships 

Pi+P2 + P3 = pl + pl + pl = i ^^ qi + q2 + q3 = ql + ql + ql = '^- (30) 

Moreover, the quantities o^a'^ in ( ^9]) are the components of a triad defining a 



Riemannian 3-dimensional geometry whose (time-independent) metric tensor is given 
by lg^,u = Saf3 oe";, oe^u, being ^e"^ the inverse matrix of oC^''. 

In this situation it is worth to note that with the results that we have obtained, 
or more precisely, from the expressions for 0, a+, cr_, and o^a^ [equations (^,^)], we 
can know the exact dependence of any quantity on the proper time r, which means that 
we have solved completely the evolution equations for the models under consideration. 
The spatial dependence of any quantity associated with these models is provided by 
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the knowledge of o^a^, oO"+, o<^-, and qT. These quantities, together with (Q) and (|29|), 
determine completely the space-time metric. It is clear that the equations for them 
come from Einstein's field equations, which in our case are equivalent to the constraints 
for the shear tensor (1Tl|), the m,om,entum, constraint: 



div(a),-|e,(e)=0, (31) 

and the equations (||) projected onto the triad {Cq,}, 

'R%^s = eA(r^,) - e,(r^,) + rs,r^, - r^,r^, - Yxs^^ = o , (32) 

where F^^ = e° ■ (Vg^^e^) are the Ricci rotation coefficients and 7^^ = e" ■ [e^, e^] are 
the commutator functions. Both quantities are associated with the triad {ea} and are 
related by 

^c^e^h = ^aeYxp + Sp.YaX + ^Xellp ■ (33) 

From (|29| ) we get the following expression for the commutator functions 

lap = {^% + {^k °^l ~ 4 °^/) [^^(^ ~ or)d^V\ - (^ - orY^VypoT 
+ ln(r - ,T')d,qx - (r - or'r\}por'\ } (r - ^ri'^-'^^-'t (r - yp~'^^-'^^_ , (34) 

where ^7^^ are the (time-independent) commutator functions associated with the triad 
{o^a\i and underlined indices do not follow the usual index summation convention. 
Now we can introduce ( p8| , |3^ into equations ( |31| , |32| ) and from them we will have to 



,e3(or + ot' - ^— [off+ - \f^o(y-\) = -i=— {{V3o(^+ - off~)o7i3 - 2o(t_ „723} ' (37) 

'J oP Vo c" 



extract time-independent equations for the quantities 060'^, oO"+, oO"-, and qT- Indeed, 
the three equations contained in (^) turn out to be equivalent to the following three 
time-independent equations 

oeiioT + ot' + -^-^) = -7=— {{V3o(r+ + oO-)^l2 + (V3o(T+ - oO-)si\z\ , (35) 

^ op \/ZoP^ ^ 

9 9 

062(0^ + or' - -z—{o(y+ + ^oO-\) = -i=— {{VSoa+ + o(^-)oll2 - 2o(7- 0723} . (36) 

SoP V^oP ^ ^ 

2 . rr .. 2 

^P 

On the other hand, the equations ( ^2]) contain six independent relationships which, 
taking into account the form of the commutators functions (0) and (p9| , |33D , are linear 
combinations (with time-independent coefficients) of terms of the form 

{r-ory^'ir-oTTY and (r - „r)"^(r - y)'"Mn"«(r - ^r) ln™^(r - ,r') , (38) 

where the exponents n^, fria are rational numbers which take only some values, for 
instance: ^3,7713 = 0, 1,2 and n^ + m^ = 2. In general, from each equation in ( P^D we 
will extract several relations for the quantities o^a^, oO'+, o<^--, and „r. To that end, 
we have carried out a careful study of the linear combinations of terms ( PB[ ) coming 
from (P^, considering the linear interdependence between the terms (^) and taking 
into account that these relations must be satisfied for all the fluid elements (labelled by 
com,oving coordinates y°') and for all the values of the proper time r. 
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To put into practice the plan just outlined, it is very convenient to consider 
two separate cases depending on whether the shear tensor is degenerate or not, or 
equivalently through (p!9D, on whether the gravito-electric tensor is degenerate or not. 
The shear tensor is degenerate if and only if one of the following quantities vanishes 

=P2-Pl, =P3-Pl, =P2-P3- 

That is to say, if and only if two Pa (or also two g^) are equal. We have analyzed both 
cases arriving at the following results: 

3.1. Algebraically general case 

In this case it is clear that all the Pa must be different. In fact, from (pOD we can 
assume, without loss of generality, that in an open domain of the space-time we have: 
Ps G (— |, 0), p2 £ (0, |), and pi G (|, 1). Studying the equations ( p2D and (|35| - |37D we 
have found that they imply: 

7^^ = , Pa = constant , qT = constant . 

Therefore, we conclude that the only possible metrics in this case are those belonging 
to the family of the Bianchi I IDMs |]12|, whose line element can be written as 

ds' = -dT^ + J2 r^'^^ir - orf^'^^'^-Kdy'^y . (39) 

a=l 

3.2. Degenerate case 

In this case two of the Pa must be equal, we can choose without loosing generality 

P2-P3 = ^^ off- = , (40) 

and from ( pO]) we can take pi = —1/3, P2 = Ps = 2/3 (gi = I, q2 = Is = 0). As a 
consequence, in this case there are fewer linearly-independent terms in ( p8D and then, 
we get also fewer restrictions. Moreover, we must take into account that the triad {oGq} 
(and also the triad {e^}) is now not completely fixed [due to (^0])], there is the freedom 
of a rotation on the 2-planes spanned by 62 and 63. Taking all these considerations 
into account, we have found that the analysis of the equations ( |5B| - P7D and ( P^ implies 



that the only possible metrics are those belonging to the Szekeres space-times [|T3| and 
such that the hypersurfaces orthogonal to the fluid velocity are flat, which were given 
14| . There are two families of solutions which in our formalism are distinguished 



m 



according to whether ei (qt) vanishes or not (in the formulation of [|T^ the key quantity 
is (3' = dr(3). 

When ei(or) = 0, the line element is given by [(y") = (x, y, z)] 



ds^ = -dr' + [(1 + Ay + Bz)t + CY r'^dx' + rs dy' + dz' , (41) 
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where A, B and C are arbitrary functions of x. And when ei(or) 7^ 0, it can be cast 
into the form 



ds= = -,,. + IlM^ 



r-„r 3 



, 2 d^ioT) 

T — „T + -- 



l2 



Sd,{\nU) 



dx^ + U^T-oT)^ dy^ + dz^ ,(42) 



where 



U = VW, W = [c 



2 I 2 

y +z 



+ 2by + 2cz + dy 



and where a, b, c, and d are any functions of x such that ad — b"^ — c^ = 1, and V and 
oT are arbitrary functions of x. 



4. Analysis and discussion 

Now, we are going to analyze the consequences of the study that we have carried out 
in the previous sections. First of all, we have determined all the IDMs in which the 
hypersurfaces orthogonal to the fluid velocity are fiat. More specifically, we have shown 
the following statement: "All the irrotational dust solutions of Einstein's equations 
with flat spatial geometry are given by the family of the Bianchi I dust solutions ^3^), 
and by the two subfamilies of the Szekeres dust models given in (^7],^^". Another 
consequence of this study is that it provides an intrinsic space-time characterization 
of these solutions. This characterization is based on the fact that within this class 
of models we can distinguish two classes according to the algebraic type of the shear 
tensor, or equivalently [through the relation (]TP|)], according to the algebraic type of 
the gravito-electric tensor, Eab- But taking into account the well-known fact that these 
solutions (^,^,^) have a vanishing gravito-magnetic tensor. Hat = 0, these two classes 
can be distinguished according to the Petrov type. Then, the intrinsic characterization 
is as follows: The Bianchi I dust models (^) are the Petrov type I IDMs with flat spatial 
geometry and the two subfamilies ( ^^2|) of the Szekeres models are the Petrov type D 
IDMs with flat spatial geometry. Both classes contain the Petrov type O models, which 
correspond to the Einstein-de Sitter space-time |jl6| . 

With regard to the Killing symmetries, in the algebraically general case the Bianchi 
I models are homogeneous cosmological models in which the 3-dimensional simply- 
transitive group of motions is Abelian. In the degenerate case the situation is more 
interesting since the models in this class belong to the Szekeres class of inhomogeneous 



solutions, and in |15| it was shown that the Szekeres solutions have no Killing vectors 
in general. We have carried out the same study for the two particular families given 
by ( p| , ^2D , and we have found that both classes of solutions have also in general no 
Killing vectors. This results reinforces the idea expressed by Collins ||^ that intrinsic 
symmetries (symmetries of submanifolds) are considerably less restrictive than space- 
time symmetries. In our case the hypersurfaces S(t) are flat, i.e., they have 6 (intrinsic) 
Killings, but we have found models without space-time Killing symmetries. 

On the other hand, as we have mention before, all these solutions (|39|j4l| , |42|) have 
a vanishing gravito-magnetic tensor, Hab = 0. Therefore, they belong to the class of 
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the silent universes |T^, IDMs with a vanishing gravito-magnetic tensor. These models 



have attracted much attention in the last years in the context of structure formation 



and their study led to conjecture [18| that all the silent universes must be either Bianchi 
I, or Szekeres, or FLRW dust models. Hence, it turns out that we have proved here this 
conjecture for the particular case in which ^R^bcd = 0. In other words we have shown 
that "All the silent universes with flat spatial geometry belong to the families of Bianchi 
I, Szekeres, and FLRW dust models". The only previous partial proof of the conjecture, 
for the vacuum case, has been given in |jl9|. Therefore, our result is the first partial 



proof of the conjecture in the non-vacuum case. 

The procedure followed in this work can be applied to other physically interesting 
classes of space-times. By one hand we can consider a more general energy-momentum 
content, specifically, irrotational perfect-fluid models with a geodesic velocity field (in 
this case the pressure depends only on the proper time r). Obviously, this includes 
the dust case with a non-zero cosmological constant. Here, it would be interesting 
to investigate whether or not these models belong to the classes of Bianchi I and 
Szafron pO[ space-times (the generalization of the Szekeres models for a perfect- 
fiuid matter content). On the other hand, we can soften the fiatness assumption by 
demanding the hypersurfaces orthogonal to the fiuid velocity to have constant curvature, 
i.e., ^Rabcd = K{T){hachbd ~ hadhbc), where K{t) denotes the Gaussian curvature. 
Since the models included in these two extensions (which can obviously carried out 
simultaneously) satisfy ( pO| , pTD , the procedure we have to follow is essentially the same. 



Another issue that would deserve more attention is the question of finding a covariant 
proof of the results presented here, specially taking into account the fact that these 
models admit an intrinsic and covariant characterization. This would consist in pushing 
forward the study made in section |^, perhaps by using computer algebra. 

Finally, we want to point out that using the covariant space-time characterization 



that we have obtained for the solutions (|39| , ^ , ^2D , it is possible to find a covariant initial 



data characterization, which consists in considering the local initial-value problem for 
IDMs and to find the initial data whose development corresponds with these models. 
However, this is beyond the scope of this paper. An analysis of the initial data for some 



classes of IDMs, including those considered here, will be presented in [^ . 
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